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^ ' Abstract 

I We discuss the phenomenology of the lightest stops in models where R-parity is 

' broken by bilinear terms. In this class of models we consider scenarios where the 

r~| ■ R-parity breaking two-body decay ti ^ r+ 6 competes with the leading three-body 

C^! decays ii W~^bxi, h H~^bxi, ^^t and ii ^ bi>il~^ {I = e^n^r). 

Q^l We demonstrate that the R-parity violating decay can be the dominant one. In 

^ ■ particular we focus on the implications for a future e+e^ Linear Collider. 

> 

ix^ ■ 1 Introduction 



The search for supersymmetry (SUSY) [|I], 0] plays an important role in the experimental 
program at the colliders LEP2 and Tevatron. It will be even more important at future 
colliders, e.g. an upgraded Tevatron, LHC, an e^e~ linear collider. Therefore many 
phenomenological studies have been carried out in recent years (see e.g. 0, 0, |^ and 
references therein). Most of them have been carried out in the context of the minimal 
supersymmetric standard model (MSSM) p, 0. However, neither gauge invariance nor 
supersymmetry requires the conservation of R-parity. Indeed, there is considerable the- 
oretical and phenomenological interest in studying possible implications of alternative 
scenarios [§] in which R-parity is broken [l^, 0, |12|. The violation of R-parity could 
arise explicitly |T3| as a residual effect of some larger unified theory or spontaneously, 



through nonzero vacuum expectation values (vev's) for scalar neutrinos j^, [IT], [T^. In re- 
alistic spontaneous R-parity breaking models there is an SU{2) ® U{1) singlet sneutrino 
vev characterizing the scale of R-parity violation [1^, [1^, [1^ which is expected to be 



the same as the effective supersymmetry breaking scale. 

There are two generic cases of spontaneous R-parity breaking models to consider. In 
the absence of any additional gauge symmetry, these models lead to the existence of a 
physical massless Nambu-Goldstone boson, called majoron (J) which is the lightest SUSY 
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particle, massless and therefore stable. If lepton number is part of the gauge symmetry 
and R-parity is spontaneously broken then there is an additional gauge boson which gets 
mass via the Higgs mechanism, and there is no physical Goldstone boson fl^. As in 
the standard case in R-parity breaking models the lightest SUSY particle (LSP) is in 
general a neutralino. However, it now decays mostly into visible states, therefore diluting 
the missing momentum signal and bringing in increased multiplicity events which arise 
mainly from three-body decays such as 

X? - //^, (1) 
where / denotes a charged fermion. The neutralino also has the invisible decay mode 

XS - 3^.. (2) 

as well as 

X? ^ i^J, (3) 



in the case the breaking of R-parity is spontaneous [T^, [T^. This last decay conserves 
R-parity since the majoron has a large R-odd singlet sneutrino component. 

Owing to the large top Yukawa coupling the stops have a quite different phenomenology 
compared to those of the first two generations of up-type squarks (see e.g. ||T8| and 
references therein). The large Yukawa coupling implies a large mixing between tL and 



tji and large couplings to the higgsino components of neutralinos and charginos. 
The large top quark mass also implies the existence of scenarios where all two-body decay 
modes of ti (e.g. ti — > t Xi, b x/j t q) are kinematically forbidden. In these scenarios higher 
order decays of ti become relevant: [^, |2l[|: ti cxi2i ~^ W'^hxi, ii H^hxi 
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ti blfui, ii buil^, where / denotes e,/i,r. In |]^ it has been shown that the 



three-body decay modes are in general much more important than the two body decay 
mode in the framework of the MSSM. Recently it has been demonstrated that not only 
LSP decays are sign of R-parity violation but that also the light stop is possible candidate 
for observing R-parity violation even if R-parity violation is small |P, 22, It has been 



demonstrated that there exists a large parameter region where the R-parity violating 
decay 

h^hr (4) 

is much more important than 

ii^cx?,2 (5) 

in scenarios where only those decay modes are possible. It is therefore natural to ask if 
there exist scenarios where the decay ti ^ 6 r is as important as the three-body decays. 
Note that in the R-parity violating models under consideration the neutral (charged) 
Higgs-bosons mix with the neutral (charged) sleptons. These states are denoted by 5"°, Pj*, 
and for the neutral scalars, pseudoscalars and charged scalars, respectively. Therefore 
in the R-parity violating case one has the following three-body decay modes: 

h ^ W+hx\ (6) 

h - Stbxl (7) 

h ^ Stbui (8) 

ti ^ 6 5°/+, (9) 

h ^ 6P°/+. (10) 
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We will demonstrate that ti br^ can indeed be the most important decay mode. In 
particular we will consider a mass range of ti, where it is difficult for the LHC to discover 



the light stop within the MSSM due to the large top background p4|. The rest of paper 



is organized in the following way: in the next section we will introduce the model. In 
Sect. ^ numerical results for stop decays are presented and their implications for LC. In 
Sect. 4 we present our conclusions. 



2 The model 



The supersymmetric Lagrangian is specified by the superpotential W given by 



b 

2 ' 



(11) 



where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU{2) indices, and e is a com- 
pletely antisymmetric 2x2 matrix, with ei2 = 1. The symbol "hat" over each letter 
indicates a superfield, with Qi, Li, Hi, and H2 being SU{2) doublets with hypercharges 
1/3, —1, —1, and 1 respectively, and U, ID, and R being SU{2) singlets with hypercharges 
— |, |, and 2 respectively. The couplings hu, h^, and He are 3x3 Yukawa matrices, and 
/i and ej are parameters with units of mass. 

Supersymmetry breaking is parametrized by the standard set of soft supersymmetry 
breaking terms 

Vsoft = MfQTQ''^ + MfU*Uj + Mfb*bj + MfLTL'; + MfR*Rj 
- [IM3A3A3 + IMA2A2 + iM'AiAi + h.c] 



(12) 



Note that, in the presence of soft supersymmetry breaking terms the bilinear terms ej 
can not be rotated away, since the rotation that eliminates it reintroduces an R-Parity 
violating trilinear term, as well as a sneutrino vacuum expectation value ||2^ . 



For our discussion it suffices to assume R-parity Violation (RPV) only in the third gen- 
eration. However we do allow for R-parity-conserving Flavour Changing Neutral Currents 
(FCNC) effects, such as the process ti — ex? involving the three generations of quarks. 
In this case we will omit the labels i,j in the soft breaking terms. In order to study 
the R-Parity violating decay mode ii —* b r it is sufficient to consider the superpotential 



W = htQ^U^H2 + hQsDsHi + KL^R^Hi - fiHiH2 + e^L^^H^ (13) 

This amounts to neglecting the effects of RPV on the two ffist families. A short discussion 
on ti ^ bl'^ in the three generation model will be given at the end of Sect. |[ 

The bilinear term in Eq. (|13|) leads to a mixing between the charginos and the r-lepton 
which in turn leads to the decay ti ^ br. The mass matrix is given by 





Mr 



M j^^gv2 



(14) 
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As in the MSSM, the chargino mass matrix is diagonahzed by two rotation matrices U 
and V 



U*McV 
























irtr 



(15) 



The hghtest eigenstate of this mass matrix must be the tau lepton {t^) and so the mass is 
constrained to be 1.7771 GeV. To obtain this the tau Yukawa couphng becomes a function 



of the parameters in the mass matrix, and the full expression is given in |29 
The stop mass matrix is given by 



kv%' + A 



UR 



(16) 



with AuL = -sig'-y'M-vl + vl) and A^,^ = Ig'^vj-vl + vl)^ 
given by = g^{v'^+vl+v^)/2. The mass eigenstates are given by ti 
and ^2 = tji cos Oj—t^ sin 9^. The sfermion mixing angle is given by 



The sum of the vf is 
= ti cos Oj+tR sin 6-j: 



cos 9j 



{Ml-mlY + {Ml 



2 \2 



sin 6j 



til 



2 

m~ 



(M| - 

til 



+ (M| )2 



;i7) 



In addition the charged Higgs bosons mix with charged sleptons and the real (imag- 
inary) parts of the sneutrino mix the scalar (pseudoscalar) Higgs bosons. The formulas 
can be found e.g. in |]29|, |3^. Their main decay modes for the mass range considered in 
this study are: 



S7 



bb, T+T , Xl^r 
bb, T^T' , Xll^r 



(18) 
(19) 
(20) 



3 Numerical results 

In this section we present our numerical results for the branching ratios of the higher order 
decays of ti. Here we consider scenarios where all two-body decays induced at tree-level 
are kinematically forbidden. Before going into detail it is useful to have some approximate 
formulas at hand [^: 

Tih^br) ^ ' ' ' (21) 

^ ^ 167r ^ ^ 

nii^cxl) ^ 10"^/iX-i (22) 

where IC/32I ~ l^s/f^l if l^sl ^ \f^\ and ^3 ^ rriw- The complete formulas are given in 
[ p^ , p^ . For the three-body decays the formulas given in can be used as a good 



approximation if the mixings induced by R-parity violation are small. The complete 
formulas for the three-body decays in the R-parity violating case will be given elsewhere 
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a) Br(t,) b) Br(t,) 
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Figure 1: Branching ratios for ti decays for m^-^ = 220 GeV, /i = 500 GeV, M = 240 GeV, 
and rriiy = 100 eV. In a) the branching ratios are shown as a function of cos for tan /? = 4, 
in b) as a function of tan/3 for cos% = 0.25. Note, that the graph - - displays the sum 

bUrSt + bT+ + bT+ Pj^. 



Input: 


tan/5 = 


= 4 


^ = 500 GeV 


M = 


240 GeV 




Mn = 


370 GeV 


Mq = 


340 GeV 


At = 


150 GeV 




Me = 


190 GeV 


Mi = 


190 GeV 


Ar = 


150 GeV 




mi, = 


220 GeV 


cos 9i 


= 0.25 


irtpo 


= 110 GeV 


Calculated 


m^o = 


120 




= 225 




= 520 




Al 




= 










m-b, = 


340 GeV 


= 


375 GeV 


cos Oi 


= 0.925 




rrigo = 


82 GeV 


nigo = 


= 128 GeV 


mgo - 


= 182 GeV 




rripo = 


110 GeV 


rripo = 


= 182 GeV 








rrio- = 


= 136 GeV 


= 

^3 


= 187 GeV 


rric- 


= 204 GeV 




mg^ = 


213 GeV 


mu, = 


m^^ = 204 GeV 







Table 1: Input parameters and resulting quantities used in Fig. 1. 



We have fixed the parameters as in |2T[ to avoid colour breaking minima: we have used 
m^^, cos9i, tan/?, and /i as input parameters in the top squark sector. For the sbottom 
(stau) sector we have fixed Mq,M^ and Ab {M^,Mi, and Ar) as input parameters. In 
addition we choose the R-parity violating parameters 63 and f 3 in such a way that the 
tau neutrino mass is fixed ([^ and references therein): 



AMM'/^ - 2(c/2M' + g'^M)^'Vuv'a cos ^ ' 



yy ivi ^ y ivi ^ ,2 . 2c ^OQ^ 



with 



sine = , + (24) 



^ = \/f^' + 4, Vd=\/vj + vl (25) 

For simplicity, we assume that the soft SUSY breaking parameters are equal for all gen- 
erations. 

In Fig. ^(a) and (b) we show the branching ratios of ti as a function of cos%. The 
parameters and physical quantities are given in Tab. |l|. In Fig. |l|(a) we show BR(ti — >■ 
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6r+), BR(ti ^ cXi), BR(ti ^ BR(ti 6e+i>e) + BR(ti ^ 6z/ee+). The 

branching ratios for decays into jli oi are practically the same as those into cl or z/g. 
We have summed up those branching ratios for the decays into sleptons that give the 
same final state, for example: 

tl ^ 6 i/e ^ 6 e+ z/e X? , ti ^ 6 e+ z>e ^ 6 e+ Xi (26) 
Note that in Fig. |I] we have also summed the decay branching ratios BR(i(:i b S'^ v^.) + 

In the above cases the assumption mi^ — rrih < m-+ implies m-+ > mj. Therefore, 
charginos can not arise as decay products of sleptons. The latter can only decay into the 
corresponding lepton plus Xi except for a small parameter region where the decay into X2 
is possible. However, this decay is negligible due to kinematics in that region. In addition 
there exists the possibility of R-parity violating decays. However, these will be small 
because the neutrinos mix mainly with higgsinos implying that the partial decay widths 
are proportional to the squared product of an R-parity violating mixing parameter and 
small Yukawa coupling. For this set of parameters BR(ti cXi) is O(10~^) independent 
of cos 6*^7 and therefore negligible. Near cosOi = —0.3 one has ti —>■ bW^ Xi dominant 
decay channel, since the ti-xt-b coupling vanishes implying that the main contribution 
for the decays into the scalars vanishes. Moreover, the width for ti ^ br^ is somewhat 
suppressed because of the cos^ 6'^ factor in Eq. (^). Note, from the figure that the 
branching ratios for the various decays into selectrons e-sneutrino is roughly a factor two 
smaller than the sum of the decays into S*^, 5]^, P^. The reason is that, for this choice 
of parameters the P2 is mainly the pseudoscalar Higgs boson of the MSSM with mass 
110 GeV. In this case the R-parity violating channel ti — > fer+Pg is comparable to the 
corresponding R-parity conserving decays. This state appears additional to the states 
which carry tau-lepton number in the MSSM limit giving rise to the observed difference. 
Note, that one has to expect additional jets from the states containing the scalars 5"°, P? 
and 5^ because they have admixtures of the original Higgs boson. In case of negative 
cosOi the decay into ti — > 6r+ is important and can be even the most important one. 
Therefore one has events with t~ bb in the final state which can be used for a full mass 
reconstruction of the light stop. 

In Fig. |l|(b) the tan (3 dependence of branching ratios is shown. For this specific choice 
of cos^f the decay ti ^ br^ is the most important one for tan/? < 15. Above this value 
the final states which contain the scalars corresponding to the lighter MSSM stau are the 
most important ones. The growth of the decay branching ratios into these states with 
tan /3 is a feature independent of cos 9^. 

The assumption that no tree-level-induced two-body decays are kinematically allowed 
implies that m-+ > — rrib- Therefore, one expects an increase of BR(ti — b Xi) if 
m^^ increases, because the decay into b Xi is dominated by the t exchange, whereas for 
the decays into scalars Xj exchange dominates. This trend is indeed observed in Fig. ^ 
where we show the branching ratios for m^^ = 350 GeV. Here we concentrate on the range 
of cos 9i where At <1 TeV to avoid possible minima in the scalar potential which break 
either color or electric charge. Notice that for the heavy stop case the decay ti ^ b Xi 
is the most important one, independently of cos 6^ and tan (3. 

Note however that also in this case R-parity violation implies a distinct signature 
compared to what is expected in the MSSM due to the decays of Xv One gets the 
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Figure 2: Branching ratios for ti decays for m^^ = 350 GeV, /i = 750 GeV, M = 380 GeV, 
and mj, = 1 keV. In a) the branching ratios are shown as a function of cos 6^ for tan /? = 4, 



in b) as a function of tan (3 for cos Of 

bUrSt + bT+ + bT+ Pj^. 



0.7. Note, that the graph - - displays the sum 



a) BR(ti ^ bn)/BR{ii br) b) BR(f^ ^ 6e)/BR(ii ^ 6r) 




Figure 3: Ratio of branching ratios: a)BR(ti 6/i)/BR(ti ^ 6r) as a function of (€2/63)^ 
and b) BR(ti ^ 6e)/BR(ti br) as a function of (ei/eg)^ for for m^-^ = 220 GeV, 
/i = 500 GeV, M = 240 GeV, and ej + 4 + = I GeV^; cos^^ = 0.05,0.1,0.2,0.40.8, 
m^^ = 100 eV. 



following high-multiplicity final states: 

h bW+ffu (27) 
^ bW+ff'l^ (28) 
^ 6W^+z/J (29) 

where / denotes a Standard Model fermion. Here the decays of the PF-boson will give 
additional leptons and jets. Therefore, one has in general additional jets and leptons 
compared to the MSSM case. 

In the event that ei,2 are of the same order of magnitude, as suggested by a solution to 



the present neutrino anomalies one has in addition the decays into be'^ and If 
one passes from the 1-generation model to the 3-generation model the situation changes 
as follows. From Eq. (|2T|) it follows that the sum of the modes r{ti bl^) in the 3- 
generation model is nearly equal to T{ti — br^) in the 1-generation model, if (e'l)^ + 
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(^2)^ + (^2)^ = ^3 where the e[ (es) are the parameters of the 3-generation (1-generation) 
model. In Fig. ^ we show the ratios of branching ratios for different ti ^ bl~^ modes 
versus the ratios of different e's squared and for different values of cosO^. In both cases 
we have fixed + + e| = 1 GeV^ and in Fig. ^a ei = £3 whereas in Fig. ^ €2 = £3. 
One can see that the dependence is nearly linear even for rather small cos 6^. This result 
depends on (e^ + €-1 + eD/fi^ and on the neutrino mass m^^, since both determine the 
mixings of the leptons with the charginos. The lines indicated in the figure come closer 
to the diagonal if {el + el + e|)//i^ increases and rriu^ decreases. 



4 Conclusions 

We have studied the phenomenology of the lightest stop in scenarios where the R-parity 
violating decay ti br^ competes with three-body decays. We have found that for 
m^^ < 250 GeV there are regions of parameter where ti — 6 r+ is an important decay 
mode if not the most important one. This implies that there exists the possibility for 
full stop mass reconstruction from r+ 6 6 final states. Moreover, in this mass range 
the discovery of the hghtest stop might not be possible at the LHC (certainly this is the 
case in the MSSM). This implies that one has to take into account the importance of 
this new decay mode when designing the stop search strategies at a future e+e~ Linear 
Collider. Spontaneously and bilinearly broken R-parity violation imply additional leptons 
and/or jets in stop cascade decays. Looking at the three generation model the decays 
into ii ^ hl^ imply the possibility of measuring e^/e^ and e^/e^ and thereby probing the 
parameters associated with the present solar and atmospheric neutrino anomalies. 
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